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Only certain galaxies are included in surveys: those bright and large enough to be detectable as 
extended sources. Because gravitational lensing can make galaxies appear both brighter and larger, 
the presence of foreground inhomogeneities can scatter galaxies across not only magnitude cuts but 
also size cuts, changing the statistical properties of the resulting catalog. Here we explore this size 
bias, and how it combines with magnification bias to affect galaxy statistics. We demonstrate that 
photometric galaxy samples from current and upcoming surveys can be even more affected by size 
bias than by magnification bias. 



I. INTRODUCTION 

Any survey using cither counts of galaxies or shapes 
of galaxies must make choices: vifhich galaxies are in and 
which are out. For example, a magnitude limited survey 
contains all galaxies brighter than a fixed flux threshold. 
Because gravitational lensing can make galaxies appear 
brighter, a background galaxy whose line of sight passes 
through an overdense region might be catapulted into the 
sample when it otherwise would not have been included 
in the catalog. As was pointed out many years ago [1-6], 
this leads to non-zero correlations between background 
and foreground objects, an effect which goes by the mon- 
icker of magnification bias [7-9]. 

A less studied phenomenon is that the increase in size 
of a background galaxy by a foreground matter overden- 
sity [10] will have a similar effect if the catalog under 
consideration includes a size cut. We refer to this addi- 
tional source of bias as size bias, and to the net effect 
of gravitational lensing on the galaxy sample as lensing 
bias. 

The purpose of this paper is to introduce size bias as 
an additional effect that needs to be considered in studies 
of galaxy statistics. Any study which correlates galaxies 
in a flux- and size-limited sample with some other ob- 
servable is potentially sensitive to size as well as mag- 
niflcation bias. In addition, size and magnification bias 
affect shear observables as well. This effect has so far not 
been studied in detail and is fleshed out in a companion 
paper [11]. 

II. LENSING BIAS 

Qualitative discussion. Figure 1 illustrates the basic 
physics behind lensing bias. The figure shows size vs. i 
band apparent magnitude for a small random subsample 
of stars and galaxies in the GOODS field [12]. Our size 
variable is r, the Full Width at Half Maximum (FWHM) 
of the galaxies in the i band. The FWHM was mea- 
sured using adaptive moments [13, 14]. All images have 
been degraded by adding noise and 0.9" seeing typical 
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FIG. 1: The distribution of galaxies (red triangles) and stars 
(blue) in the GOODS field as a function of size and mag- 
nitude, after downgrading images to an assumed 0.9" see- 
ing, typical of near-future large imaging surveys such as DES 
and Pan-STARRS. The thick vertical and horizontal lines de- 
note our fiducial magnitude and size cuts (ias < 24 and 
r > 1.2 arcsec, respectively). The hatched regions contain 
galaxies included in the survey due to lensing by a foreground 
mass: the vertical band corresponds to the traditional magni- 
fication bias effect, while the horizontal region contains galax- 
ies affected by the size bias introduced in this work. "Stars" 
at lower left with abnormally small sizes are shot noise. 



of the upcoming Dark Energy Survey [DES, 15]. Also 
shown in the figure are horizontal and vertical blue lines 
corresponding to a size and magnitude cut of r > 1.2" 
and iAB < 24, respectively. The magnitude cut is nec- 
essary to obtain high S/N images in surveys with DES- 
like depth, while the size cut is meant to be comfortably 
larger than the 0.9" seeing of the experiment, which is 
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necessary in order to be able to robustly separate galax- 
ies from stars, and to obtain good shape measurements 
of the galaxies^. 

The diagonal dotted lines in Figure 1 are contours 
of constant surface brightness. Since lensing preserves 
surface brightness, a foreground matter overdensity will 
move all sources along these diagonal lines upwards and 
to the right, increasing the number of objects in the 
galaxy sample. Conversely, an underdense region along 
the line of sight can act in the opposite way. In the weak 
lensing limit, the flux magnification factor A is given by 
A ~ 1-|-2k, where k is the convergence (see below). Then, 
the magnitude and size of a source will be (de-)magnified 
according to: 
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The narrow horizontal and vertical hatched bands shown 
in Figure 1 contain galaxies which would pass our galaxy 
selection cuts if they were to be lensed by matter along 
the line of sight. For illustration purposes, we chose a 
large value of the lensing convergence, k = 0.1. The well 
known magnification bias effect corresponds to the ver- 
tical hatched region. In addition, however, we see that 
there is a comparable population of galaxies in the hor- 
izontal hatched region which is brought into our galaxy 
sample not because the galaxies are made brighter, but 
because the galaxies become more extended. 

Quantitative treatment. We wish to determine how the 
galaxy density field nobs(^) observed in an experiment re- 
lates to the intrinsic (i.e. un-lensed) galaxy density field 
uq^O). Our discussion follows closely the presentation 
found in the appendix of Hui ct al. [16]. Let /, r, and 
6 denote the observed flux, apparent size, and position 
of a galaxy, respectively, and fg, Tg, and 6g denote the 
corresponding intrinsic quantities. Further, let ^{f,r,9) 
and ^g{fg,rg,9g) denote the observed and intrinsic dis- 
tribution of galaxies in flux, size, and position. Because 
the total number of galaxies is preserved, one has that 



$(/, r, e) df dr d^e = ^g{fg,rg, Og) dfg drg d^ 



(2) 



In the linear regime, the intrinsic and observed galaxy 
properties are related via^ 



f^Afg, r=VA 
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^ For faint objects in the downgraded data, photon noise leads to 
large scatter in the estimated rpwHM- Note that the star-galaxy 
separation is technically done based on the resolution parameter 
-R [13], not rpwHlvl, and the relation between the two has some 
scatter. 

^ Linear distortions of galaxies adequately describe gravitational 
lensing everywhere except where galaxies fall near the caustic of 
a foreground lens. Also, while the impact of gravitational lensing 
on the source radius r may depend on the precise definition of r, 
we expect the relation r = A^I'^Vg to hold in the vast majority 
of cases to a good approximation. Other scalings can be trivially 
included following the discussion presented in this work. 



where (56' is the gravitational lensing deflection angle, and 
A = \dBldBg\ is the corresponding amplification. We 
work in the weak lensing regime throughout, so that A = 
1 -I- 2k. Note that the surface brightness of a source S oc 
//r^ is unaffected by A. 

Let now e( f, r) be the galaxy selection function of a 
survey, which we assume depends only on a source's mag- 
nitude and size. Given this selection function, the intrin- 
sic and observed galaxy number densities are: 
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M^g) = j '^fa j '^^9^if9^rg)<S>g{fg,rg,9g). (s) 

Using equation (3), we can write /, r, and 6 in the above 
expression in terms of the corresponding intrinsic quanti- 
ties. Linearizing, and using equation (2) to replace $ by 
$3, we find that the relation between the observed and 
intrinsic galaxy densities is given by 

nobs(el = no{0)[l + {2(3 f + /?. - 2)k], 
where the parameters /?/ and (3r arc defined via: 
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While strictly speaking (3r are defined in terms of 
the intrinsic galaxy distribution $g, we have replaced 
this with the observed galaxy distribution $, which is 
the only one accessible to observations. Any differences 
between the two are of order k averaged over the survey 
area, and hence correspond to a higher order correction 
in our perturbative approach. 

If we further assume that the selection function is given 
by a simple magnitude and size cut, the function e(f,r) 
takes the form e(/,r) = e(/ - /min)6(r - r^in) where 
Q{x) is a step function. In this case, the expressions for 
f3f and Pr simplify to 
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Note that in this case, one has /3f ~ {h/2)d\gn/dm = 
5s/2. If Pr = 0, one obtains riobs = "^o[l + (5s — 2)k], 
which is the standard expression for magnification bias. 
For Pr Q, size bias acts in a way that is analogous to 
magnification bias, except that it is sensitive to the slope 
of the galaxy size distribution rather than that of the 
magnitude distribution. The parameter which captures 
the full extent of lensing bias therefore is 



q=2Pf + Pr- 2. 



(10) 



To first order in k, lensing bias affects the observed num- 
ber density via 



(11) 
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FIG. 2: Magnitude (ias, top panel) and apparent size (Igr, 
bottom panel) distribution of galaxies in the GOODS field. In 
each case, the solid red line shows the full-resolution GOODS 
data, while the blue dashed line shows the distribution for 
data degraded to the expected DES instrument, typical noise, 
and 0.9" seeing. 



FIG. 3: Flux- and size count slopes, 2/3/ and measured 
from GOODS data degraded to DES seeing and noise, as a 
function of the magnitude cut for a fixed size cut (r > 1.2 arc- 
sec; top panel), and as a function of the size cut for a fixed 
magnitude cut (iab < 24 mag; bottom panel). Also shown 
is the total q from Eq. (11) (solid line), and the q obtained 
when neglecting size bias (long-dashed line). 



The next section presents estimates of q for upcoming 
galaxy surveys. 

III. FORECAST FOR UPCOMING SURVEYS 

Figure 2 shows the galaxy density of the GOODS field 
as a function of both magnitude and size. The GOODS 
data is degraded to 0.9" seeing and typical noise, and 
the galaxies are error weighted using adaptive moments 
[13; 14]. In order for galaxy ellipticities to be reliably 
measured, one requires both {ab ^ 24 and r > 1.2", 
which we adopt as our fiducial magnitude and size cuts. 

In order to measure the q parameter for this data 
set, we determine /?/, /3r from the downgraded GOODS 
galaxy sample according to Eq. (9), for a range of cut 
values. Note that, consistent with our perturbative ap- 
proach, we have neglected the intersection between the 
two hatched regions in Fig. 1 in this calculation, which 
would further increase the value of q but is formally of 
order k^. 

While we take into account realistic weights assigned 
to each galaxy from the measurement errors, the results 
presented here are mainly for illustrative purposes: sharp 
size and fiux cuts are rarely made in practice. Neverthe- 
less, in many cases some form of apparent size cut is ap- 
plied to galaxy catalogs. For instance, the main observ- 
able of weak lensing studies is galaxy ellipticity, which 
can only be reliably measured if the galaxies are signif- 
icantly more extended than the PSF of the instrument. 



Consequently, weak lensing studies usually introduce a 
size and magnitude dependent error weighting function 
e(m, r) which operates as a selection function. In general, 
this selection will be a smooth function of m, r, rather 
than a sharp cut. In that case, the resulting q will be 
a suitable average over the sharp-cut q values presented 
here. 

Figure 3 shows the parameters f3 / and Pr as a function 
of both the magnitude and size cuts, while holding the 
other cut fixed. The corresponding value of q is shown 
(solid line) as well as the q obtained when neglecting size 
bias (long-dashed line). It is clear from the figure that 
size bias can have a significant impact on q for size cuts 
larger than the instrumental PSF, eventually dominat- 
ing the lensing effect for large cuts. The fact that size 
cuts smaller than the PSF have no impact on the total 
lensing bias is not surprising: since all galaxies have an 
apparent size at least as large as the PSF, there is ef- 
fectively no size cut if Tcut is smaller than the PSF. We 
conclude that the sample of source galaxies in the DES 
weak lensing survey will be affected by size bias at a level 
that is comparable to or larger than the corresponding 
magnification bias. Interestingly, we find that for typical 
cut values, the q parameter including size bias is positive, 
while it is negative when no size cut is applied. It might 
be possible to confirm this prediction with existing data 
from galaxy surveys. As shown in the companion paper 
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[11] and briefly in Sec. TV, lensing bias will turn out to be 
significant when interpreting the results from upcoming 
cosmic shear experiments such as DES. 

We have also investigated whether size bias could en- 
hance the lensing signature in the Luminous Red Galaxy 
(LRG) sample in the Sloan Digital Sky Survey, and in 
particular in the correlation function of these galaxies. 
The sample docs include an effective size cut in order to 
separate stars from galaxies. However, since these galax- 
ies are quite bright and large, most of them lie far from 
the thresholds. Therefore, size bias plays an insignificant 
role in the interpretation of the LRG correlation function. 
We estimated that this will also hold for higher-redshift 
LRG samples derived from the SDSS. For a different rea- 
son, size bias is expected to be unimportant for quasars 
(QSO): since these are point-like objects, lensing docs 
not affect their size noticeably, so that QSO are again 
only affected by magnification bias. 

IV. EFFECT OF LENSING BIAS ON SHEAR 
MEASUREMENTS 

As an application of our discussion of lensing bias, we 
briefly present the effect on shear observablcs which is 
the subject of [11]. Since the cosmic shear field 7 can be 
measured only where there are background galaxies, 7 is 
preferentially sampled in regions with a high density of 
background galaxies. However, a part of the fluctuations 
in the galaxy density come from lensing bias which is due 
to the same lensing field as the shear 7 itself. Because 
of this, the actual measured shear field is schematically 
given by (see [11] for details): 

ioUe) = [i + qm]^{0). (12) 

The leading lensing bias correction to the shear power 
spectrum is formally identical to the reduced shear 
correction [17, 18]. However, the amplitude of the 
correction is multiplied by a factor 1 -I- g « 2 — 3, from 



our estimates in Sec. HL Thus, the correction to the 
shear power spectrum from lensing bias is at the 5% 
level for I ^ 1000, and growing for higher £ [11]. 

V. SUMMARY AND CONCLUSIONS 

Gravitational lensing magnifies galaxies, making them 
not only brighter, but also larger in appearance. If galaxy 
selection criteria include any form of size measure, it fol- 
lows that gravitational lensing can scatter galaxies in and 
out of the sample across the applied size cut. This ef- 
fect, which we have dubbed size bias, is completely anal- 
ogous to magnification bias, and is important whenever 
a size cut is applied that removes some non-negligible 
fraction of galaxies from the sample. Size bias will of- 
ten be just as important as magnification bias, so many 
previous examples subject to the latter might profitably 
be re-examined for traces of size bias. Our calculations 
in §HI suggest that upcoming large galaxy surveys will 
be similarly afflicted by both types of lensing bias. In a 
companion paper [11], we evaluate the impact of lensing 
bias on shear observables, and demonstrate that this ef- 
fect is significant at the 5% level for I ^ 1000 in case of 
the shear power spectrum. We estimate that for DES, 
neglecting lensing and size bias can lead to biases in the 
Dark Energy parameters estimated from cosmic shear at 
the 2—3(7 level. Similar conclusions are expected to apply 
to other upcoming surveys such as Pan-STARRS [19]. 
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